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We give a detailed description of the entanglement purifi- 
cation protocol which generates maximally entangled states 
with high efficiencies from realistic Gaussian continuous vari- 
able entangled states. The physical implementation of this 
protocol is extensively analyzed using high finesse cavities 
and cavity enhanced cross Kerr nonlinearities. In particular, 
we take into account many imperfections in the experimental 
scheme and calculate their influences. Quantitative require- 
ments are given for the relevant experimental parameters. 
PACS numbers: 03.67.Hk, 42.50.-p, 03.65. Bz 

I. INTRODUCTION 

Quantum entanglement plays an essential role in many 
interesting quantum information protocols, such as in 
quantum key distribution and quantum teleportation [Q . 
To faithfully realize these protocols, first we need to gen- 
erate a maximally entangled state. In reality, however, 
due to loss and decoherence, normally we can only gen- 
erate partially entangled states between distant sides . 
Entanglement purification is further needed which distills 
a maximally entangled state from several pairs of par- 
tially entangled states using local quantum operations 
and classical communications [^1). For qubit systems, 
efficient entanglement purification protocols have been 
found [HJ^. Recently, quantum information protocols 
have been extended from qubit systems to continuous 
variable systems, such as continuous variable teleporta- 
tion 1^,1^ , continuous variable computation |^ and error 
correction , continuous variable cryptography , and 
also the notions of continuous variable inseparability ||l^] 
and bound entanglement have been investigated. For 
physical implementation, Gaussian continuous variable 
entangled states (i.e., states whose Wigner functions are 
Gaussians) can be generated experimentally by transmit- 
ting two- mode squeezed light, and this kind of entangle- 
ment has been demonstrated in the recent experiment of 
continuous variable teleportation Obviously, it is 

useful to consider purification of continuous variable en- 
tanglement, that is, to generate a desired more entangled 
state from some realistic continuous entangled states. We 
have recently proposed an efficient continuous variable 
entanglement purification protocol p^ . In this paper, 
we present the mathematical details of this purification 
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protocol together with new results on its physical imple- 
mentation. In particular, we take into account many im- 
portant imperfections in a realistic experimental setup, 
and calculate their influence on the purification scheme. 
Quantitative requirements are given for the relevant ex- 
perimental parameters. These calculations make neces- 
sary preparations for a real experiment. We also show 
how to generate Gaussian continuous entangled states 
between two distant high finesse cavities, which is the 
first step for the physical implementation of the purifica- 
tion protocol. 

It should be noted that with direct extensions of the 
purification protocols for qubit systems, it is possible to 
increase entanglement for a special class of less realis- 
tic continuous entangled states . Unfortunately, with 
these direct extensions no entanglement increase has been 
found till now for realistic Gaussian continuous entangled 
states. In [|l^ a protocol to increase the entanglement 
for the special case of pure two-mode squeezed states 
has been proposed, which is based on conditional photon 
subtraction. For its practical realization, the efficiency, 
however, seems to be an issue. In contrast, the purifi- 
cation scheme discussed in this paper has the following 
favorable properties: (i) For pure states it reaches the 
maximal allowed efficiency in the asymptotic limit (when 
the number of pairs of modes goes to infinity) ; (ii) It can 
be readily extended to distill maximally entangled states 
from a relevant class of mixed Gaussian states which re- 
sult from losses in the light transmission; (iii) An experi- 
mental scheme is possible for physical implementation of 
the purification protocol using high finesse cavities and 
cross Kerr nonlinearities. 

The paper is arranged as follows: In section 2 we 
show how to generate a Gaussian continuous entangled 
state between two distant cavities from the broadband 
squeezed light provided by a nondegenerate optical para- 
metric amplifier (NOPA). Light transmission loss is taken 
into account. In sections 3 and 4 we give a detailed de- 
scription of the purification protocol. Section 3 shows 
how to generate a maximally entangled state from pure 
two-mode squeezed states based on a local quantum non- 
demolition (QND) measurement of the total photon num- 
ber, and section 4 extends the purification protocol to 
include the mixed Gaussian continuous states which are 
evolved from the pure two-mode squeezed states due to 
the unavoidable light transmission loss. In section 5, we 
describe a cavity scheme to realize the local QND mea- 
surement of the total photon number, and deduce condi- 
tions for the QND measurement. Then, in section 6, we 
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extensively discuss many imperfections for a real exper- 
iment on QND measurements, and deduce quantitative 
requirements for the relevant experimental parameters. 
Last, we summarize the results, and give some typical 
parameter estimations. 



II. GENERATION OF CONTINUOUS 
ENTANGLED STATES BETWEEN TWO 
DISTANT CAVITIES 

Our source of entangled light field is taken to be a 
NOPA operating below threshold The light fields 

may be nondegenerate in polarization or in frequency. 
The two NOPA cavity modes and cb are assumed to 
have the same output coupling rate Kc- The dynamic in 
the NOPA cavity is described by the Langevin equations 
(in the rotating frame) (l8| 
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where e is the pumping rate with |e| < Kc/2 (below 
threshold), and CiA and CiB are vacuum inputs. The 
NOPA outputs CoA and Cqb are given respectively by 
Coa = Cia + ^/K^Ca {a = A,B). The two outputs, per- 
haps after a long distance propagation, are incident on 
distant high finesse cavities A and B. The cavities A and 
B are assumed to have the same damping rate k with 
K <ti Kc- The schematic setup is shown by Fig. 1. 
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FIG. 1. Schematic setup for generating Gaussian continu- 
ous entangled states between two distant cavities. 

Under the condition k <C Kc: the dynamics in the 
NOPA cavity is much faster than those in the cavities 
A and B, so we can assume a steady state for the NOPA 
outputs. The steady NOPA outputs are described by 
squeezed white noise operators with the following corre- 
lations 

{coA (t) CoB (i')> = M5 {t - t') , 

(cL {t) Coo {t')) ^NS{t- t') , {a = A,B), (2) 
{coo. (t) cl^ it')) ^{N + l)S{t- t') , {a^A,B), 

where N and M, satisfying M = y^N (N + 1), 
are determined by the NOPA coupling and pumping 

rates through N = k^/ ^-^^ — |e|^^ and M = 

kkc(# + |er)/(4-Nf . 

To get the steady state of the cavities A and B, we 
note that their inputs qia and a^s are respectively the 



NOPA outputs CoA and Cob with neglect of the losses 
during light propagation. The Langevin equations for 
the cavity modes a a and Ub have the form 



\/Kaia, {a^A,B), 



with the following solution 
aa (t) = Qa (O)e"'^* - 



e-*(*-*')a,„(i')dt'. (3) 



When hd is considerably larger than 1, from Eqs. (2) and 
(3), it follows that 



[aACiB ) = (N+l), 
(aia^) ^N, {a = A,B), 
(a^ai) - (Af + 1), (a = A,B), 



(4) 



On the other hand, we know that two modes driven by a 
white noise are in Gaussian states at any time. A Gaus- 
sian state with the correlations (4) is necessarily a pure 
two-mode squeezed state. So the steady state of the cav- 
ity modes a a and as is 
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squeezing 



(5) 



opera- 



and the squeezing 



parameter r is determined by coth(r) = \/N + 1. 

Next we include some important sources of noise in the 
state generation process. The noise includes the losses in 
the NOPA cavity and the light transmission loss from 
the NOPA cavity to the cavities A and B. With a small 
loss rate rjo <C Kc for the modes ca and cb in the NOPA 
cavity, the Langevin equation (1) is replaced by 



CA = ec'r 
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(6) 



where ViA and ViB are standard vacuum white noise, and 
the NOPA outputs are still given by Coa = + ^/kZcci 
[a ~ A,B). On the other hand, the transmission loss of 
light can be described by 



a^a = CoaVe^^ + VaVT~^^^, (a = A, B) , (7) 

where t is the transmission time, r]A and tjb are respec- 
tively the transmission loss rates for the outputs Cqa and 
CoBi and VA and vb are standard vacuum white noise. 
From Eqs. (6) and (7), it follows that the inputs for the 
cavities A and B have the following correlations 



{a^A (t) a^B {t')) = ^N' (N' + l)e 



(aL (i) it')) = N'e-<-6 [t - t') , {a^A,B), 
(a,„ {t) (t')) = (n'c-'^'^- + l)5{t~ t') , [a ^ A, B) 
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where the total loss rates 77^ = J?a + 7 In (1 + tjq/kc) = 
Va + Vq/ (i^ct) (q: = a, and the parameter N' = 

|e|'(«c + r?o)'/(^^^^T^-kl')' ~ N. The steady 
state of the two cavity modes ua and as is thus a Gaus- 
sian state with the non-zero correlations given by 



{uAaB) = \/F(7VTT)e 



(ala„> =Ne-<^, {a = A,B), (8) 
- fNe-i'o'-' + 1) , ia = A,B). 



The Gaussian state is completely determined by these 
correlations. The Gaussian state (8) can be equivalently 
described as the solution at time t = t oi the following 
master equation 

, f tit If 
P =Va [aApa'^ - -a'^aAiP - ^pa\aA 



If If 
-a^j^asp- -pa'^ttB 



(9) 



with the initial state p (0) = 1^')^^ (*|, where 1^')^^ 
is defined by Eq. (5). This equivalence simplifies the 
physical picture in section 4, where we will use the master 
equation (9) to describe the state generation noise. 



III. ENTANGLEMENT CONCENTRATION OF 
PURE TWO-MODE SQUEEZED STATES 

In the above, we have shown how to generate continu- 
ous partially entangled states between two distant cavi- 
ties. In the case of no noise in the state generation pro- 
cess, the cavities are in a pure two-mode squeezed state. 
In this section, we will show how to concentrate continu- 
ous variable entanglement, that is, starting from several 
pairs of continuous entangled states, we want to get a 
state with more entanglement through only local oper- 
ations. The section is divided into two parts. The first 
part describes the purification protocol for two entangled 
pairs, and the second part extends the protocol to include 
multiple pairs. 



A. Concentration of two entangled pairs 

Assume now we have two cavities AijA^ and Bi, B2 on 
each side. Each pair of cavities Ai, Bi (i = 1,2) are pre- 
pared in the state (5), which is now denoted by \'^)a b-- 
expressed in the number basis, has the form 



the entanglement is uniquely quantified by the von Neu- 
mann entropy of the reduced density operator of its one- 
component. The entanglement of the state (10) is thus 
expressed as 

E {\^)a b ) ~ cosh^ r log (cosh^ r) — sinh^ r log (sinh^ r) . 

(11) 

The joint state of the two entangled pairs Ai,Bi and 
A2,B2 is simply the product 

I*)aiBiA2S2 = '^^iSi l^ac)^^s^ O SA2B2 (r) Ivac)^^^^ 



= (i-A')E^'ViT7bU^,B,B.. (12) 
3=0 

where |j)AiA2SiB2 defined as 
1 ^ 

\j) A-^A^BiB^ = y/T^ E ~ "^^1^2 l"''^' ~ '^)biS2 • 



n=0 



(13) 



We now perform a local QND measurement of the total 
photon number of the two cavities Ai,A2. There have 
been several proposals for doing QND measurements of 
the photon number, and in section 5, we will describe a 
cavity scheme for realizing the QND measurement of the 
total photon number of two local cavities. Here we simply 
assume this type of measurement can be done. After 
the QND measurement of the total number n^i + nA2 , 
the state \'^)a,b,A2B2 collapsed into \j)A,A2B,B2 ^i*^ 
probability 



(14) 



The state \j) A1A2B1B2 ^ maximally entangled state be- 
tween the two parties Ai,A2 and Bi,B2 in a (j -|- 1) x 
{j + l)-dimensional Hilbert space, and its entanglement 
is 



^^(li)AiA2BiS2) =l0g(i + l)' 



(15) 



2B1B2) -^(l*)AiBi)' l-^- j > 

(sinh(r))"°'''''> ~ 1' S^* a two-party state with more 

entanglement. The quantity F, = '^(^•^]^i'^2^i^2) defines 

the entanglement increase ratio. Fig. 2 shows the prob- 
ability of success versus entanglement increase ratio for 
some typical values of the squeezing parameter. 



(10) 



n=0 



where A = tanh(r). Equation (10) is just the Schmidt 
decomposition of the state I'S')^.^.. For a pure state, 
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FIG. 2. The purification success probability versus entan- 
glement increase ratio for two pairs. Dotted line for the 
squeezing parameter r — 0.5, dashed line for r — 1.0, and 
solid line for r = 1.5. 

An interesting feature of this entanglement purification 
protocol is that with any measurement outcome j ^ 0, 
we always get a useful maximally entangled state in some 
finite Hilbert space, though the entanglement of the out- 
come state \j} A1A2B1B2 "^^^^ necessarily exceed that 
of the original state \'^)AiBi J small. The state 
\j)AiA2BiB2 involves two pairs of cavities. If one wants 
to transfer the entanglement to a single pair of cavity 
modes, one can make a phase measurement of the cavity 
mode ^2- There have been some proposals for doing a 
phase measurement p^ , pO[ | . A phase measurement of the 
mode A2 with the measurement outcome (j) will convert 
the state \ j) A1A0B1B2 following maximally entan- 

gled state of a single pair of cavity modes 



1 



A1A2 



E 



^i{j~n)4> I 



(16) 



ji=0 



B. Concentration of multiple entangled pairs 

The above protocol can be extended straightforwardly 
to simultaneously concentrate entanglement of multiple 
cavity-pairs. Simultaneous concentration of multiple en- 
tangled pairs is much more effective that the entangle- 
ment concentration two by two. Assume that we have m 
cavity-pairs Ai,Bi, ^2,52, • • • and Am,Bm- Each pair 
of cavities Ai,Bi is prepared in the state (10). The joint 
state of the m entangled pairs can be expressed as 



1^) 



A2B2 



I*; 



A„^B„ 



^l-A^)^^AV/j™)|j)(^.Ba' (17) 

3=0 



where {AiBi} is abbreviation of Ai,Bi, ^2,^2, 
Ajyi, Bjn, and the normalized state \j)(A B } defined 



(A,B,} 



p(m) 



E 



^1 7 ^2 : ' ' ' 1 

il + ?'2 H h im = j 



(18) 



The function /j™-* in Eqs. (17) and (18) is given by 
P(m) _ (i + m-1)! _ f j + m-l 



J 1 



j\ (m — 1)! 



m — 1 



(19) 



To concentrate the entanglement, we perform a QND 
measurement of the total photon number ua^ + nA2 + 
• • ■+nA-^- This measurement projects the state \^)(^a b } 
onto a two-party maximally entangled state \i){^A B } 
with probability 



pM = (l_A2)"A2-'/ 



(m) 



(20) 



The entanglement of the outcome state \ j)(^A B } given 
by 



E 



(A^B^}) =log(/: 



Am) 



(21) 



Similarly, F^- = £; /£; (1*)^^^^ defines the 

entanglement increase ratio, and if Tj > 1, we get a more 
entangled state. For four pairs, the probability of success 
versus entanglement increase ratio is shown in Fig. 3. 
There appears a peak in the probability curve for some 
entanglement increase ratio between 2 and 3. 
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FIG. 3. The purification success probability versus entan- 
glement increase ratio for the number of pairs m — 4. Dot- 
ted line for the squeezing parameter r — 0.5, dashed line for 
r — 1.0, and solid line for r = 1.5. 
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To measure how efficient the scheme is, we define the 
entanglement transfer efficiency T with the expression 



T = 



j=0 ^ 



(22) 



It is the ratio of the average entanglement after con- 
centration measurement to the initial total entanglement 
contained in the m pairs. Obviously, T < 1 should always 
hold. With the squeezing parameter r = 0.5, 1.0 or 1.5, 
the entanglement transfer efficiency versus the number of 
pairs m is shown in Fig. 4. 
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FIG. 4. The entanglement transfer efficiency versus the 
number of pairs m in simultaneous concentration. Dotted 
line for r — 0.5, dashed line for r — 1.0, and solid line for 
r = 1.5. 

From the figure, we see that the entanglement transfer 
efficiency is near to 1 for a large number of pairs. In fact, 
it can be proven that if m goes to infinity, with unit prob- 
ability we would get a maximally entangled state with 
entanglement mE (l^')^.^.) . To show this, we calculate 

the mean value and the variance of the distribution pj™'' , 
and find 



J = 



{NT 



(1-A2)^ 



(23) 



The results show that if m tends to infinity, 
and the distribution p^™^ tends to a 5-like function. 
Furthermore, around the mean value j, the entanglement 
of the resulting state is 

so the entanglement transfer efficiency tends to unity. 
This proves that the purification method described above 



is optimal in the asymptotic limit (m —foo), analogous 
to the purification protocol presented in |Q] for the qubit 
case. For any finite number of entangled pairs, this pu- 
rification protocol is more efficient than that in Q| , since 
it takes advantage of the special relations between the 
coefficients in the two-mode squeezed state. 



IV. ENTANGLEMENT PURIFICATION OF 
MIXED GAUSSIAN CONTINUOUS ENTANGLED 
STATES 

The assumption of noise-free preparation of partially 
continuous entangled states is not realistic. If we include 
the unavoidable light transmission loss and the NOPA 
cavity loss in the state generation process, in section 2 we 
have shown that we would get a mixed Gaussian contin- 
uous entangled state between two distant cavities. The 
state is described by the solution at the transmission time 
T of the master equation (9), with the ideal two-mode 
squeezed state (10) at the beginning. If we want to es- 
tablish m entangled cavity-pairs Ai, _Bi, A2, i?2, ■ • ■ and 
Am,Bm, Eq. (9) can be extended directly to the follow- 
ing form 



P= -i H^sp - pH[ 



1=1 

(25) 



where p is the density operator of the whole m entan- 
gled pairs with p{0) = ^^'^ effective 
Hamiltonian 



E 



Vb t 



(26) 



In Eqs. (25) and (26), we assumed that the total loss 
rates ry^ and rj^ are the same for the m entangled pairs, 
but rj'j^ and 77^ may be different from each other. In this 
section, we will show how to distill entanglement from 
the kind of realistic continuous entangled states described 
by the solution of the master equation (25). There are 
two practical circumstances in which our entanglement 
purification protocol can be extended straightforwardly 
to generate maximally entangled states from the mixed 
Gaussian entangled states. We describe these two cir- 
cumstances one by one. 



A. Case of small state preparation noise 

Though the state preparation noise is unavoidable, in 
many cases it is reasonable to assume that it is quite 
small. We take ry^r and rj'gT as small factors, and solve 
the master equation (25) perturbatively to the first order 
of these small factors. It is convenient to use the quantum 
trajectory language to explain the perturbative solution. 
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In this language, to the first order of ry^r and 7]^t, the 
final normalized state of the m entangled pairs is either 
(no jumps) 



oo 




with probability 

or (a jump occurred) 



(28) 



(1 - A^)"" X'^^ f':'"''\-{lA+^B)^i. 



(32) 



We also note that the projection operators P^'' P^'' can- 
not eliminate the state obtained from the initial state 
s } by ^ quantum jump on both sides A and B. 
The total probability for this kind of quantum jumps to 
occur is proportional to m?rfiri'j^r]'gT^ . So the condition 
'■^'febSnddll^stalte preparation noise in fact requires 



p(m) 



TO n [rjAT + 77o/kc) (r/sr + ?7o/kc) < 1- 



(33) 



If the light transmission loss is the dominant noise, Eq. 
(33) reduces to m^n^ijArjBT'^ ^ 1- 



B. Case of asymmetric state preparation noise 



^VaTCLa, , {a ^ A, B and 



with probability 



>0 - 



(29) 



(30) 



where 



mat. 



I* 



'{A B } ~ sinh^ (r) is the 



mean photon number for a single mode. 

Similar to the pure state case, we also use QND mea- 
surements of the total photon number to distill entan- 
glement from the mixed continuous state described by 
Eqs. (27)-(30). The difference is that now we perform 
QND measurements on both sides A and B. The mea- 
surement results are denoted by Ja and Jb , respectively. 
We then compare Ja and Jb through classical communi- 
cation, and keep the outcome state if and only if Ja = Jb- 
It is easy to show that the final state is a maximally en- 
tangled state in a finite dimensional Hilbert space. Let 
P^ and Pg denote the projections onto the eigenspace 

of the corresponding total number operator ^ a^.a^i 

2 = 1 

m 

and ^ a^ GBi with eigenvalue j, respectively. From 
Eqs. (27) and (29), it follows 



^A ^B 



pi'}) pU) 

^A ^B 



= \3) 



{A,B,} 



0, {a = A,B and i — 1, 2, • • • , m) 



So if jA — Jb, the outcome state is maximally entan- 
gled with entanglement log ^/j™''^ The components 

(29) in the mixed density operator, which are not max- 
imally entangled, are discarded through confirmation of 
the two-side measurement outcomes. Compared with the 
pure state case, the probability to get the entangled state 



i Jnlt^^. abow) purification protocol, we need classical 
communication (CC) to confirm that the measurement 
outcomes of the two sides are the same, and during this 
CC, we implicitly assume that the storage noise for the 
cavity modes is negligible. In fact, that the storage noise 
during CC is much smaller than the transmission noise is 
a common assumption made in all the entanglement pu- 
rification schemes which need the help of repeated CCs 
[^,||. If we also make this assumption for continuous 
variable systems, there exists a simple purification pro- 
tocol to generate maximally entangled states. We put 
the NOPA setup on the A side. After creation of ideal 
squeezed vacuum lights, we directly couple one output 
light of the NOPA to the cavity on side A without noisy 
propagation; and the other output of the NOPA is sent to 
the remote side B, through a long distance noisy trans- 
mission. This configuration of the setup is equivalent 
to setting the transmission loss rate rjA ^ so that 
v'a ~ Vo/ (kct). Note that the NOPA cavity loss rate 
770 is normally much smaller than the output coupling 
rate Kc, so the total loss rate 77^ can be much smaller 
than r]'g in this case. The purification protocol now is 
exactly the same as that described in the previous case. 
We note that the component of the mixed density opera- 
tor which is kept the projection Pj['Pg' should subject 
to the same times of quantum jumps on each side A and 
B. We want this component is a maximally entangled 
state. This requires that the total probability for A and 
B to be subjected to the same nonzero number of quan- 
tum jumps should be very small. From Eq. (30), this 
total probability is always smaller than ynnry^r, no mat- 
(^•'^now large the transmission loss rjBT is. So the working 
condition of the protocol in the asymmetric transmission 
noise case is 



is now decreased to 



mnr]o/Kc ^ 1. 



(34) 



The transmission loss rjBT can be above one. The prob- 
ability of success for obtaining the maximally entangled 
state \ j)^A B } ^1^° given by Eq. (32). 
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Before concluding this section, we remark that for con- 
tinuous variable systems, the information carrier is nor- 
mally light, and the assumption of storage with a very 
small loss rate is typically unrealistic. It is interesting 
to note that recently there have been proposals to store 
light in internal states of an atomic ensemble If 
this turns out to be possible, the storage time for light 
can be greatly increased. Anyway, as was pointed out 
in jl^], this purification method is in fact not essentially 
hampered by the difficulty to store light, since there is a 
simple posterior confirmation method to circumvent the 
storage problem. Note that the purpose to distill maxi- 
mally entangled states is to directly apply them in some 
quantum communication protocol, such as in quantum 
cryptography or in quantum teleportation. So we can 
modify the above purification protocol by the following 
procedure: right after the cavity A attains its steady 
state, we make a QND measurement of the total exci- 
tation number on side A and get a measurement result 
jA- Then we do not store the outcome state on side 
A, but immediately use it (e.g., perform the correspond- 
ing measurement as required by a quantum cryptogra- 
phy protocol). During this process, the modes Bi are 
being sent to the distant side B, and when they arrive, 
we make another QND measurement of the total excita- 
tion number of the modes Bi and get a outcome Jb ■ The 
resulting state on side B can be directly used (for quan- 
tum cryptography for instance) if Ja = js , and discarded 
otherwise. By this method, we formally get maximally 
entangled states through posterior confirmation, and at 
the same time we need not store the modes on both sides. 



V. QND MEASUREMENTS OF THE TOTAL 
PHOTON NUMBER OF SEVERAL CAVITIES 

The QND measurement of the total photon number 
plays a critical role in our entanglement purification pro- 
tocol. There have been some proposals for making a 
QND measurement of the photon number in a single cav- 
ity [p3|-p5|, such as letting some atoms pass through the 
cavity, and measuring the internal or external degrees 
of freedom of the atoms In this section, we pro- 

pose a purely optical scheme for making a QND mea- 
surement of the total photon number contained in sev- 
eral cavities. The different optical modes interact with 
each other through cross phase modulation induced by a 
Kerr medium, and we use cavities to enhance this kind of 
interaction. As an illustrative example, in the following 
we will show how to measure the total photon number of 
two cavities. Extension of this scheme to include several 
cavities is straightforward. 

The schematic setup is depicted in Fig. 5. We want to 
make a QND measurement of the total photon number 
rii + 712 contained in the good cavities I and II, whose 
damping rate k is assumed to be very small. The cavi- 
ties I and II, each with a Kerr type medium inside, are 
put respectively in a bigger ring cavity. The two ring 



cavities are assumed to damp at the same rate 7, and 
7 ^ K. A strong coherent light bn is incident on the first 
ring cavity, whose output boi is directed to the second 
ring cavity. The output 602 of the second ring cavity is 
continuously observed through homodyne detection, and 
we will show that under some realistic conditions, this 
detection gives a QND measurement of the total photon 
number operator ni + n2 — a\ai + 0302. 




FIG. 5. Schematic experimental setup to measure the total 
photon number ni -I- n2 contained in the cavities I and II. 



The measurement model depicted in Fig. 5 is an ex- 
ample of a cascaded quantum system [p^.The incident 
light bii can be expressed as 

6^1+5^7, (35) 

where 17^/7 is a constant driving field, and b^ is the stan- 
dard vacuum white noise, satisfying 

{b?At)Kiit')) =0, 

{Ki{t)b'^At')) ^6{t-t'). (36) 

The Hamiltonian for the Kerr medium is assumed to be 

H, = hxn,b\b,, (i = l,2), (37) 

where 61 and 62 are the annihilation operators for the 
ring cavity modes, and x is the cross-phase modulation 
coefficient. The self-phase modulation effects will be dis- 
cussed in the next section and shown to be negligible 
under some realistic conditions. In the rotating frame, 
the Langevin equations describing the dynamics in the 
two ring cavities have the form 

bi = -ixnibi - -bi - V7&-i - gj, 

7 

62 = -ixnibi - - \/lba (38) 

The boundary conditions for the two ring cavities are 
described by 



7 



bi2 = boi = 6-1 + g,/y + ^/yh, 
bo2 = bi2 + V762. 



(39) 



Assume 7 ^ x{i^i) : (* = lj2), and we take adiabatic 
elimination, i.e., let bi—b2— in Eq. (38), obtaining 



bi 



-2 {gi+Vlb'a) 

7 

2 (g7 + V7^^i) 
7 



1 - 



1 - 



7 

4ixni 
7 



2^X'^ts 
7 



(40) 



Substituting the above result into Eq. (39), the final 
output field bo2 is expressed as 



'Jo2 



' v/7 



(m + 712) + 6-1 +5^7- 



(41) 



Now wc mcasiirc the X-component of the quadrature 
phase amplitudes of the output field 602 through a ho- 
modyne detection. The phase of the driving field g is set 
according to g = i\g\. Suppose T is the measuring time. 
What we really get is the integrated photon current over 
time T, which, divided by T, corresponds to the following 
measuring operator 



1 



1 



T 



4V2 



72 



bo2 it) + &I2 (*) 



9\X 



dt 



.(6) 



%/7 



(ni+n2) + --^X^ , 



(42) 

where x!^^ = '^{^T + b^x^, and 6t, satisfying 



bT,b\ 



= 1, is defined by 



(43) 



From Eq. (36), it follows that the defined mode bx is in 
a vacuum state. So the first term of the right hand side 
of Eq. (42) represents the signal which is proportional 
to ni + 712, and the second term represents the contribu- 
tion of the vacuum noise. The distinguishability of this 
measurement is given by 



Sn= 

U 6n < 1, i.e., if the measuring time 

7 



64|5rx'' 



(44) 



(45) 



we perform an effective measurement of the total number 
operator ni+n2- During the measuring time T, the loss 
of the two cavities I and II should be negligible, which 
requires 



K{ni)T< 1, (i = l,2) 



(46) 



Under this condition, rii + 712 is approximately a con- 
served observable, and we realize a QND measurement 
of the total photon number operator. The measurement 
projects the field in the cavities I and II to one of the 
eigenstates of rii +n2- Eqs. (45) and (46), combined to- 
gether, determine the suitable choice for the measuring 
time. 



VI. INFLUENCE OF IMPERFECTIONS IN THE 
QND MEASUREMENT 

We have shown how to perform a QND measurement 
of the total photon number. The scheme described above 
works under ideal conditions. For a real experiment, 
there are always many imperfections which should be 
considered. For example, the phase of the driving field 
may be unstable, and has a small variance; the damp- 
ing rates and the cross phase modulation coeSicients for 
different ring cavities may not be exactly the same; the 
Kerr media and the mirrors may absorb some light; self- 
phase modulation effects caused by the Kerr media may 
have some influence on the resulting state; there may be 
some loss of light from the first ring cavity to the second 
ring cavity; the efficiency of the detector is not unity. Of 
course, to realize a QND measurement of the total pho- 
ton number, all the imperfections must be small. But 
the important question is how small these imperfections 
should be. In this section, we will deduce quantitative re- 
quirements for all the imperfections listed above. These 
calculations may be helpful for a future real experiment. 
We will consider these imperfections one by one. 



A. Phase instability of the driving field 

Assume that the phase of the driving field g^^ has a 
small variance S, i.e., g is expressed as g = i \g\ e*'^. Then, 
Eq. (42) is replaced by 



(b) 



-y/2\g\5^, 

(47) 



The last term of Eq. (47) represents the noise due to 
the phase instability of the driving field. It should be 
negligible compared with the signal, which requires 



^ < 



4x 
7 ' 



(48) 



On the other hand, we know that the squared phase 
variance increases linearly with time, i.e. = 6tt, 
where St is the increasing rate. The measuring time T is 
bounded from below by Eq. (45), so the increasing rate 
of the phase instability of the driving field is required to 
satisfy 
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St< 



1024|g|'x^ 



(49) 



Eq. (49) suggests it is easier to meet the requirement 
imposed by the phase instabihty with a strong driving 
field and a large cross phase modulation coeflacient. 



and Cii and Ci2 arc the standard vacuum inputs. The 
boundary conditions for the ring cavities are still de- 
scribed by Eq. (39). The leaked (or absorbed) light fields 
Col and Co2 are expressed as 



+ ^/Kba, (a = 1,2). 



(55) 



B. Imbalance between the ring cavities 

In the previous section, we assumed that the damping 
rates and the cross phase modulation coefficients arc ex- 
actly the same for the two ring cavities. This may be 
impossible in a real experiment. Here we calculate the 
largest allowed imbalance between the two ring cavities. 
The damping rates and the cross phase modulation co- 
efiicients for the ring cavities are denoted by 71, 72 and 
Xii X2, respectively. The Langevin equations (38) and 
the boundary conditions (39) are replaced respectively 
by the following equations 



lib'ii -911, 



71 

bi = -ixmih - 

72 

62 = -«X2"'2^2 - Y^2 - Vl2bi2 



bi2 = boi = &■! + + Vti^I' 

bo2 = bi2 + 1/7^62- 

The final measured observable is expressed as 



(50) 



(51) 



The leakage (or absorption) of light may have two types 

of effects: First, it may destroy the balance between the 
two ring cavities; and second, the leaked light (55) may 
carry some information about rii (or 1x2). Any informa- 
tion about ni (or 712) will destroy the superposition of 
the different eigenstates of ni (or 712), and thus lead to 
decoherence of the eigenstate of rii + 112 (Note that a 
eigenstate of ni -|- 712 is normally a superposition of the 
different eigenstates of n\ (or 112)). So we require that 
the information about t?,i (or 77.2) carried by the leaked 
light should be completely masked by the vacuum noise. 
This is equivalent to require that the decoherence of the 
eigenstate of 77,1 -I- 712 caused by the light leakage is neg- 
ligible. To consider the first effect of the light leakage, 
we calculate the measured observable Xt, and find it has 
the form 



Xn 



4^2 



9\X 



1 



(6) 



V7 



(m + 712) + -^Xj. 



+ 



Ay/2 



9\X 



/3| 



PI 



(56) 



The last term of Eq. (56) should be negligible compared 
with the signal, which requires 



Xj 



Xi 



(m + 712) + -^^P + 4^/2 \g\ VtT 



, X2 

71 — 
72 

(52) 



The last term of Eq. (52) represents the noise due to 
the unbalance between the ring cavities, which should be 
negligible compared with the signal, yielding 



X271 



1 



X172 



< 



1 



(«2)- 



(53) 



C. Absorption and leakage of the light 

Light absorption by mirrors and Kerr media and light 
leakage through other mirrors of the ring cavities can be 
described by the same Langevin equation, which has the 
form 



7 / — 

61 = -ixnibi - -bi - yfyb'ii - 57 - yfei - y/PiCa, 

7 02 I — 

62 = -ixn2b2 - - ^/lbi2 - —b2 - v/^2Ci2, (54) 

where /3i and /?2 are the light leakage (or absorption) 
rates of the first and second ring cavities, respectively. 



/3i 



Xi , 
— "2, 

71. 



|/3|-/3?|< 



7 



(57) 



To consider the decoherence effect of the light leakage, we 
define a similar measuring operator X^^ for the leaked 
fight (55) 



X. 



^ 8V2 \g\xVK {a -1) 



Coa (t) + cL (t)] dt 



7 



(ni + ri2) + ^#"^ (58) 



AV2\g\xVK 



7 



n„, (a = 1,2), 



where Xj^"\ 



similar to Xj^^ defined below Eq. (42), 



are standard vacuum noise terms. The last term of Eq. 
(58) bears some information about n^, which should be 
completely masked by the vacuum noise term to make 
the decoherence effect negligible. This condition requires 



4\/2 |5f| X (no) 



7 



//3a < 



1 



2T 



(59) 



On the other hand, the measuring time T is bounded 
from below by Eq. (45), which, combined with Eq. (59), 
yields the following requirement for the leakage rates 
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f3a < 



1,2) 



(60) 



Obviously, this is a much stronger requirement than that 
given by Eq. (57). 

We should mention that there is another kind of ab- 
sorption by the Kerr medium, the absorption rate of 
which is proportional to the cavity photon number tIq. 
This kind of absorption, usually termed two-photon ab- 
sorption, cannot be described by Eq. (54). To incor- 
porate the two-photon absorption, we add an imagi- 
nary part to the cross phase modulation coefficient 
i.e., X is replaced by x -I- ixi, where Xi describes the 
two-photon absorption rate. The two-photon absorption 
should be negligible compared with the cross Kerr inter- 
action, which requires Xi < "(Try' ("^ = 1, 2) • 



D. Self-phase modulation effects 

Normally, a Kerr medium also induces self-phase mod- 
ulation effects. However, by a suitable choice of the reso- 
nance condition for the Kerr medium, the self-phase mod- 
ulation effects can be made much smaller than the cross- 
phase modulation then the self phase modulation 
interaction is basically negligible. Here, for complete- 
ness, we still calculate the influence of self-phase modu- 
lations. In fact, self-phase modulation of the ring cav- 
ity modes have no influence on the QND measurement. 
This modulation adds a term like —iXsb\hihi [i = 1, 2) in 
the Langevin equation (38), where Xs denotes the self- 
phase modulation coefficient for the ring cavity modes. 
We know that the ring cavity modes hi and 62 are in 
steady states under adiabatic elimination, and to a good 

approximation h\bi can be replaced by ib\bi^ — 4|(/|^. 

So the term —ixsb\hibi simply induces a constant phase 
shift for the output field 602, and it can be easily com- 
pensated by choosing the initial phase of the driving field 
9- 

Self phase modulation of the cavity modes ai and 02 
plays a more subtle role. First, it obviously has no in- 
fluence on the QND measurement of ni + n2, but it in- 
fluences the resulting state after the QND measurement. 
In the purification scheme for two entangled pairs (de- 
scribed in section III. A), if there is no self-phase modu- 
lation, the state after the QND measurement is given 
by Eq. (13); and if the self-phase modulation of the 
modes ai and 02 is considered, the modulation Hamil- 
tonian hx'snf {i = 1,2), in which x's is the corresponding 
self phase modulation coefficient, will bring the resulting 
state into 



1 

UIA1A2B1B2 ~ /Tin ^ ' '-^ 'MA2 

(61) 



where t is the interaction time for the self phase modula- 
tion. It is important to note that the state (61) is still a 
maximally entangled state with entanglement log {j + 1). 
In this sense, self-phase modulation effects have no influ- 
ence on the entanglement purification, though the result- 
ing state is changed. 



E. Imperfect coupling from the first ring cavity to 
the second ring cavity 

If the coupling between the two ring cavities is not 
perfect, the relation bi2 = 60I is not valid any more, and 
should be replaced by 



do = \/J^di + VT~~^5oi, 



(62) 



where di is the standard vacuum white noise, and do rep- 
resents the leaked light in the imperfect coupling. The 
quantity /i describes the coupling efficiency. This kind of 
imperfection is very similar to the light leakage (or ab- 
sorption) described in subsection VI. C. The difference is 
that the imperfect coupling (62) does not cause any un- 
balance between the two ring cavities. The only restric- 
tion is that the decoherence effect induced by it should 
be negligible, which requires 



^J.> 1- 



1 



(63) 



Eq. (63) suggests that loss of light from the first to the 
second ring cavity should be very small. 



F. Detector inefficiency 



The detector efficiency of course cannot attain 1. For 



a detector with efficiency v, the real measured field 6J, 



o2 



has the following relation with the output of the second 
ring cavity 



K2 = V^bo2 + VT^ 



(64) 



where is the standard vacuum white noise. This im- 
perfection is similar to the imperfect coupling considered 
in the previous subsection. But now the leaked light de- 
pends only on the operator sum ni + n2, and carries no 
information about the single cavity photon number m, so 
it does not induce any decoherence. The only role played 
by the detector inefficiency is that it decreases the sig- 
nal by a factor y/i/, so Eq. (45) on the restriction of the 
measuring time is now replaced by 



T> 



7 



6Av\g\ x^ 



(65) 



Obviously, the detector inefficiency has no important in- 
fluence on this QND measurement scheme. 
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VII. SUMMARY AND DISCUSSION 

In summary, we have given a detailed description of 
the purification protocol which generates maximally en- 
tangled states in a finite dimensional Hilbert space from 
two-mode squeezed states or from realistic Gaussian con- 
tinuous entangled states. The nonlocal Gaussian contin- 
uous entangled states are generated by feeding two dis- 
tant cavities with the outputs of the NOPA. The purifi- 
cation operation is based on a local QND measurement 
of the total photon number contained in several cavi- 
ties. We have extensively analyzed a cavity scheme to 
do this QND measurement, and have deduced its work- 
ing condition. Furthermore, we have discussed many im- 
perfections existing in a real experiment, and deduced 
quantitative requirements for the relevant experimental 
parameters. In Tab. 1, we summarize the working condi- 
tions for the collective QND measurement, including the 
requirements for many types of imperfections. 



Measuring time 



Phase instability 



Cavity imbalance 



Absorption (leakage) rate 



Coupling efRciency 



Detector efficiency 



<T < 



jC ^ 4y c - 1024|orx 

< ^ or Ot < '3 ^ 



^271 _ 1 < 1 

X1T2 ("2) 



/3c. <j^, (a = 1,2) 



At > 1 



1 



v > 



TABLE I. List of requirements for the QND measurement 
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To realize the QND measurement, basically we need 
high finesse optical cavities and strong cross Kerr in- 
teraction media. A good example for the strong cross 
Kerr interaction is provided by the resonantly enhanced 
Kerr nonlinearity, which has been predicted theoretically 
p6| , ^ and demonstrated in recent experiments ||2^. In 
those works, the Kerr medium is a low density cold 
trapped atomic gas, whose relevant energy level struc- 
ture is represented by the four-state diagram shown in 
Fig. 6 with |1) being the ground state. The ring cav- 
ity mode bi with frequency is assumed to be resonant 
with the |1) — |3) transition, and the cavity mode 
with frequency uja {uJa is quite different from ujb) is cou- 
pled to the |2) |4) transition, but with a large detuning 
A42 . A nonperturbative classical coupling field with fre- 
quency LOc resonant with the |2) — > |3) transition creates 
an electromagnetically induced transparency (EIT) for 
the cavity fields ai and bi . 




enough for performing QND measurements of the pho- 
ton number. For example, if the mean photon number 
(ni) — (712) — sinh^ (r) ^ 1.4 with the squeezing param- 
eter r ^ 1.0, we choose the decay rates k/2tt ~ 4MHz 
and 7/27r ~ lOOMHz (these values for decay rates are 
obtainable in current experiments ), and let g ^ 50 (for 
a cavity with cross area S ~ 0.5 x 10~^cm^, g ^ 50 cor- 
responds to a coherent driving light with intensity about 
lOmWcm^^). With the above parameters, all the re- 
quirements listed in Tab. 1 can be satisfied if we choose 
the measuring time T ~ 8ns. Note that the lig ht speed 
can be much reduced in the EIT medium p8| , so it is 
possible to get a reduced cavity decay rate n with the 
same finesse mirrors, and then more favorable parame- 
ters can be given for the QND measurement. Note also 
that a large Kerr nonlinearity based on EIT can also be 
obtained in other systems, such as trapping a single atom 
in a high finesse cavity So the example discussed 

here is not the unique choice. 
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GG acknowledges support by the Friedrich-Naumann- 
Stiftung. 



FIG. 6. Level structure of the atoms. 

In this configuration, the one-photon absorption of the 
medium is eliminated due to quantum interference, and 
the cross Kerr nonlinearity is only limited by the two- 
photon absorption (the self Kerr nonlinearity is negligi- 
ble provided that |cJq — cj^l ^ ^42)- After adiabatically 
eliminating all the atomic levels, the cross phase modu- 
lation coefficient is given by pq] 



3 \gi3f |g24|^ 
X ~ A 'T-atonn 



(66) 



where (724 and 513 are the coupling coefficients between 
the atoms and the cavity modes and bi, respectively, 
flc denotes the Rabi frequency of the coupling field, and 
riatom is the total atom number contained in the cavity. 
The two-photon absorption rate Xi is connected with x 
by the relation Xi/x = 742 /A42, where 2742 is the spon- 
taneous emission rate from level |4) to level |2). To justify 

the adiabatic elimination, one requires that 



Igial "aton 



< 



1 



30[ . As an estimation, if one takes ^sial^^ta 



0.2, 

524/27r ~ lOMHz, 742/27r ~ 30MHz, and A42 ~ IO742, 
the coefficient x is about x/'^'^ ~ 0.2MHz, and the two- 
photon absorption rate Xi ~ 0.1%. This value of the 
cross phase modulation coefficient x is not large enough 
to realize a single-photon turnstile device |2fl], but it is 
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